Test of mutually unbiased bases for six-dimensional photonic quantum systems 
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In quantum information, complementarity of quantum mechanical observables plays a key role. 
If a system resides in an eigenstate of an observable, the probability distribution for the values 
of a complementary observable is flat. The eigenstates of these two observables form a pair of 
mutually unbiased bases (MUBs). More generally, a set of MUBs consists of bases that are all 
pairwise unbiased. Except for specific dimensions of the Hilbert space, the maximal sets of MUBs 
are unknown in general. Even for a dimension as low as six, the identification of a maximal set of 
MUBs remains an open problem, although there is strong numerical evidence that no more than 
three simultaneous MUBs do exist. Here, by exploiting a newly developed holographic technique, 
we implement and test different sets of three MUBs for a single photon six-dimensional quantum 
state (a "qusix"), encoded either in a hybrid polarization-orbital angular momentum or a pure orbital 
angular momentum Hilbert space. A close agreement is observed between theory and experiments. 
Our results can find applications in state tomography, quantitative wave-particle duality, quantum 
key distribution and tests on compleinentarity and logical indeterminacy. 

PACS numbers: 



INTRODUCTION 

The "complementarity" of different observables of a 
same physical system is one of the basic features of the 
quantum world [T] . Besides its fundamental relationship 
with the concept of wave-particle duality, complementar- 
ity today plays a key role in areas ranging from quantum 
state reconstruction [2] to quantum information appli- 
cations such as quantum key distribution |2|. Mathe- 
matically, complementary observables are described by 
noncommuting Hermitian operators whose sets of eigen- 
states form different bases in the Hilbert space that are 
said to be "mutually unbiased" (MUBs) . This expression 
refers to the fact that the overlap (or inner product) of 
any pair of states belonging to different bases is the same 

HI- 

In quantum cryptography, complementary observables 
and the associated MUBs are the core of all protocols 
proposed for secure quantum key distribution, starting 
from the famous BB84 protocol |S] and its extension to 
three qubit bases [6,. The "no cloning theorem" implies 
that Alice and Bob can always recognize a possible eaves- 
dropper attack by detecting the associated disturbance 
introduced in the system. In particular, the adoption 
of MUBs for encoding the information is known to maxi- 
mize this disturbance allowing one to recognize the attack 
most effectively [3j. 

In high-dimensional systems, complementary observ- 
ables and the corresponding MUBs have been exploited 
to enhance the security in quantum cryptograpy [7_ , per- 
form fundamental tests of quantum mechanics, such as 
quantum contextuality [51I10|. explore logical indetermi- 



nacy |11| . and several other tasks in quantum informa- 
tion. For instance, new quantum key distribution pro- 
tocols were conceived in which a larger error rate can 
be tolerated while preserving security |12l [T^ . More- 
over a different protocol extending Ekert91 [14] by us- 
ing entangled qutrits has been experimentally realized 
I IS ] . In quantum state tomography. MUBs play a crucial 
role because they correspond to the optimal choice of the 
measurements to be performed in order to obtain a full 
reconstruction of the density matrix. 

Given a Hilbert space of dimension d, an important 
problem is to find the maximum number of MUBs that 
can be defined simultaneously. Although for spaces of 
prime-power dimensions there exist several methods to 
find a maximal set of d -I- 1 MUBs |T5] , this problem re- 
mains in general hitherto unsolved |17l 118] . Dimension 
six, in particular, has been widely investigated in the last 
few years |19ff22] because it is the lowest one for which 
the problem is still open. Nevertheless, several theorems 
imposing restrictions on the properties of the maximal set 
of MUBs have been proved |23l ^^ and strong numeri- 
cal evidence suggests that no more than three mutually 
unbiased bases actually exist in dimension six |19 } [22 | f25]. 

Different experimental approaches have been recently 
adopted to implement complete sets of MUBs for state 
reconstruction in photonic systems. For example, the po- 
larization of a photon pair was used to define MUBs in di- 
mension four [53], the orbital angular momentum (0AM) 
of single photons has been used to address Hilbert spaces 
with d — 2,3,4,5 |27| . and multiple propagation modes 
were combined to reach dimensions d = 7, 8 j28| . Hy- 
brid methods combining polarization and 0AM were also 



used to define and manipulate photonic ququarts {d — 4) 
|29L [501 • However, since in d = 6 no complete set of 
MUBs is known, this case has not been investigated hith- 
erto for state tomography and, to our knowledge, even 
the minimal set of three MUBs has never been demon- 
strated in an experimental framework. In the present 
work, we shall focus our attention on this important six- 
dimensional (6D) case. 

In this paper we demonstrate the generation of MUBs 
in 6D by exploiting two different approaches. A 6D 
Hilbert space can be always decomposed in the ten- 
sor product of a two-dimensional (2D) and a three- 
dimensional (3D) space. Hence, a possible route to im- 
plementing 6D quantum states (qusix) is by combining 
the 2D and 3D spaces related to two different degrees of 
freedom of the photon. Following this route, in our first 
experiment we prepare and analyze all states of three 
MUBs in a 6D Hilbert space obtained by combining the 
2D space of polarization and a given 3D subspace of the 
0AM. The 0AM degree of freedom is related to the pho- 
ton's transverse-mode spatial structure, and in principle 
it allows one to define a single-photon Hilbert space of ar- 
bitrarily large dimensionality ^3Tk34| . In our second ex- 
periment, we exploit this feature to prepare and test three 
MUBs in a 6D space defined only in the 0AM of the pho- 
ton. The polarization qubit is easily manipulated with 
standard optical elements, while arbitrary 0AM states 
can be obtained by computer-generated holography us- 
ing a spatial light modulator (SLM). In the present work 
we developed a novel method to determine the hologram 
kinoform in order to obtain a sufficiently high fidelity in 
the state generation. Indeed, MUBs are very sensitive to 
the generation fidelity and relatively small imperfections 
are immediately visible in the MUBs cross overlaps. A 
high fidelity is also crucial to exploit MUBs for quantum 
cryptography. 



MUTUALLY UNBIASED BASES 

Let A and B be two operators in a d-dimensional 
Hilbert space, with orthonormal eigenbases {|ai)} and 
{|6i)} respectively. Eigenstates of these observables are 
said to be mutually unbiased |2l |4] if 
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Such operators are also called mutually complementary, 
or maximally noncommutative, since given any eigen- 
state of one, the outcome resulting from a measure- 
ment of the other is completely undetermined. In a d- 
dimensional Hilbert space a pair of MUBs can always be 
found. Indeed, let {|a,i)} be the computational basis. 



A discrete Fourier transform can be then used to define 
the following dual basis, which is mutually unbiased to 
the previous one: 
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where cj^ ~ exp (i27r/(i) , and the non-italic i denotes the 
imaginary unit (not to be confused with the index i). 
The pair of operators associated to these bases, often 
named Z and X reminiscent of the Pauli operators, pro- 
vides an algebraic complete set of observables that fully 
parametrizes the physical degree of freedom described by 
the Hilbert space: all other operators acting on this space 
are product of powers of Z and X |1] . 

An open issue concerns the maximal number of MUBs 
that can be found in a d-dimensional space; in the spe- 
cific case when d is equal to a prime number or to a 
prime power, a maximal set of d+1 MUBs does exist |3]. 
This set is also "complete", in the sense that by projec- 
tive measurements over its states (d— l)(d+ 1) = d^ — 1 
independent real parameters can be obtained, which are 
exactly the number of parameters needed for full den- 
sity matrix reconstruction |2 • A complete set of MUBs 
can be found using several methods, i.e., the Galois Field, 
the Heisenberg-Weyl group, Hadamard matrices, etc. (for 
a review see |31 US]). However, in the general case of 
composite dimensions that are not prime powers such as 
d = 6, 10, 12, ..., all these methods fail [35]. On the base 
of extensive numerical simulations, it has been conjec- 
tured that complete sets of MUBs do not exist in this 
case [25|, although such conjecture hitherto has not been 
rigorously proven. A minimum number of MUBs that is 
known to exist in such cases is given by p'^ + 1, where 
p'' is the lowest factor in the prime decomposition of the 
number d [3S]- For instance, in the d = 6 case, three 
MUBs can be easily constructed, but no evidence for the 
existence of a fourth basis that is unbiased with the first 
three has ever been found. 

The Hilbert space H^ of a 6D system can be always 
factorized in the direct product of a 2D and a 3D space. 



I.e. 



W ^H' 



^W 



(4) 



Both these Hilbert subspaces possess complete sets of 
MUBs, {|wf)} and ^p?/f, containing respectively 
three and four bases. Although the states of these bases 
can be combined to form twelve different separable bases 
for the space T-L^ , which can be used for a complete to- 
mography of the qusix state, only sets with three MUBs 
can be constructed. A possible choice is given by the 
following three bases: 



{|a,)} = {|0),|l),...,|d-l)}. 
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where i G {1,2} and j G {1,2,3}; an explicit matricial 
expression of states that we consider is reported in the 
Methods section. It can be immediately seen that any 
other basis obtained introducing the fourth one {["•I)} 
of the space "Ha would not be mutually unbiased with the 
others, since it is missing a different basis in 7i2- This 
set of 18 product states cannot be extended by any other 
vector in He, even if entangled states are considered ^23j : 
moreover if a complete MUB set in d = 6 existed, then 
only one among the seven bases therein could be com- 
posed of product states, while all others must be entan- 
gled Pl]. 



OAM ENCODING AND THE HOLOGRAPHIC 
TECHNIQUE 

Let I to) denote the eigenstates of the photon OAM 
with eigenvalue mh, where to is any integer. These 
states define the entire infinite-dimensional Hilbert space 
of OAM. We note that a full specification of the trans- 
verse optical modes would actually require assigning also 
a radial number, as in the case of Laguerre-Gauss modes. 
However, here and in the following we will omit this ra- 
dial number and the specification of a given radial profile 
will be understood for each value of \m\. 

We now define a 3D subspace O in OAM as that 
spanned by the three eigenvectors {|-l-l) , |0) , | — 1)}. 
These states can be taken to define the logical basis 
Oi of a photonic qutrit in the Hilbert space O. A sec- 
ond basis in O that is mutually unbiased with the log- 
ical one can be obtained as the Fourier-transform one 
©2 = {|ai) , |a2) , las)} defined as in Eq. ([3|. Explicitly, 
we have |ai) = (|-l) + |0) + |l))/y3, jaa) = (|-l)+w|0)- 
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|+l))/y3, lag) = (1-1) + uj^ |0) + uj |+l))/\/3, where 
= a;3 = exp (i27r/3) and we have used the identity 
Lj" — uj. The other two bases of a maximal set of MUBs 
in O are then defined as follows: O3 = {|/3i) , j/^s) , l/^a)} 
with |/3i) - (1-1) +c^|0) +c^|l))/V3,|/?2) = (1-1) + 
6.2 |0) + |+l))/y3, 1/33) = (1-1) + |0) + ^2 |+l))/,/3 and 
Oi = {|7i),l72),l73)} with |7i) = (1-1) +Lo^\0) + 
^^|1))/V3,|72) = (1-1) +a.|0) + |+1))/^,|73) = 
{\-l) + \0)+co\+l))/V3. 

As shown in Eq. (Isl , in order to construct three MUBs 
in the 6D hybrid space, we will only need the first three 
bases 01,02,03- The intensity and phase profiles of 
the nine OAM states belonging to these three MUBs are 
shown in Fig. [l] 

Let us now discuss the experimental method we 
adopted for the generation (and detection) of these nine 
states, and of all other OAM superposition states in this 
work. Arbitrary optical field transverse modes can be ob- 
tained by diffraction of an input Gaussian TEMqo mode 
on a SLM programmed for displaying a prescribed "kino- 
form", that is the pattern determining the phase retar- 
dation experienced by the input wave in diffraction. The 




Figure 1: MUBs for hybrid photonic qusix encoding: 

Representation of quantum states with dimension d — 6 ob- 
tained from the direct product of a three-dimensional sub- 
space O of OAM and the two-dimensional space tt of polar- 
ization. The three main boxes correspond to the three MUBs. 
On the left side, the intensity and phase distributions of each 
OAM spatial mode and the corresponding generating kino- 
form are shown. On the right side the polarization states are 
illustrated graphically by showing the optical electric field ori- 
entation in space at a given time. 



main problem is that the SLM is a phase-only optical 
element, while to obtain arbitrary OAM modes we need 
to be able to tailor both the phase and the amplitude 
transverse profiles of the outgoing field. This can be ob- 
tained by modulating both the shape and contrast of the 
kinoform fringes. To determine the kinoform, we initially 
tried some of the most commonly used methods |37j . but 
found that they often give rise to a non-negligible "cross- 
talk", i.e., nonzero overlaps between different states of 
the same basis, and to significant unbalances in the over- 
laps of each state of a given basis with the states of other 
bases. In other words, the generation fidelity of these 
methods was not good enough for our purposes. For 
this reason, we developed an holographic method that 
is specifically optimized in the fidelity. 

Let us first assume that the input field is a plane wave. 
Our goal is to obtain in the first order of diffraction a 
prescribed optical field ^e**, A and $ being the optical 
field normalized amplitude and phase, respectively. A 
straightforward calculation [39J shows that such optical 
field is obtained in the far field if the kinoform phase 



modulation has the following expression: 



M ^ Mod 
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where A is the grating period that fixes the diffraction 
angle, X — (l + sinc~^(^)/7r), in which sinc~^ stands 
for inverse of sinc(x) = sin(a;)/a; function in the domain 
[— TT, 0], and Mod is the modulo function that gives the 
remainder after division of the first argument by the sec- 
ond. The inverse of sine function was evaluated numeri- 
cally by the Newton method, with an accuracy of seven 
digits. 

By this method we calculated the kinoforms needed to 
generate the nine 0AM states of the first three MUBs 
of the 0AM qutrit. The resulting hologram patterns are 
shown in Fig. ^ It can be seen that these kinoforms 
include only an azimuthal dependence, since the 0AM 
state definition ignores the radial coordinate. This im- 
plies that the same kinoforms can also be used with a 
Gaussian input beam instead of a plane wave and only 
the radial profile of the diffracted wave will be affected, 
while the 0AM state will remain the same. Moreover, 
we do not need to finely adjust the input beam waist of 
the Gaussian beam. 

We note that the holograms defined by Eq. ^ gen- 
erate ideally exact modes in the far field, so that the 
expected overlap between states belonging to the same 
basis vanishes identically and that between states be- 
longing to different MUBs is 1/3 in the qutrit space (and 
hence it will be 1/6 in the qusix space, after combining 
with polarization). As mentioned, this is not the case for 
other commonly used holographic methods. For exam- 
ple, numerical simulations based on the method reported 
in Ref. [38, yielded mean state fidehties of 100%, 88.5% 
and 84.5% for the three 0AM qutrit bases. This corre- 
sponds to 6% and 7.7% of mean cross-talk between states 
belonging to the last two bases. Moreover, the overlap 
between states of different MUBs is found to vary be- 
tween 21% and 48%, depending on the state pair. These 
fidelity problems are absent in our method. More details 
about the performances of the holographic method used 
in this work will be reported elsewhere [39] . 



HYBRID QUSIX ENCODING AND 
CHARACTERIZATION 



Our first experimental implementation of qusix pho- 
tonic states has been achieved by combining the 2D space 
TT of polarization and the 3D subspace O of 0AM in sin- 
gle photons. The logical basis / of quantum states in 
dimension six has been hence implemented as follows: 

/ = {\H,-l) , \H,0) , \H, +1) , \V, -1) , \V,0) , \V, +1)} 

(7) 



where \H) , \V) denote horizontal and vertical linear po- 
larizations, as shown in Fig.^ which have been combined 
with the three eigenstates of 0AM forming the basis Oi. 
Following Eq. ^, a second basis //, unbiased with the 
first, is obtained by combining the diagonal/antidiagonal 
polarization states {|^) = T^^^'' + 1^))' 1^) = 
-^i\H) - \V))} with the 0AM states of basis O2. The 
third basis of the set of MUBs was finally obtained 
by combining the circular polarization states {\L) — 
j^{\H) + z \V)),\R) = ^(|i/) - z \V))} with the 0AM 
third basis O3. 

In order to experimentally generate these hybrid qusix 
states we employed the setup shown in Fig. ??. Single 
photons emitted at a 7 kHz rate via spontaneous para- 
metric down conversion (SPDC) in a beta-barium-borate 
(BBO) crystal [40, are collected by a single-mode fiber 
(SM) to filter out all the spatial modes but the Gaussian 
mode TEMoo, i.e., 0AM state |0) (0AM qutrit initial- 
ization). A set of waveplates (C) compensates the po- 
larization after the transmission through the fiber. The 
photons are then sent through a polarizing beam splitter 
(PBS) (polarization qubit initialization) and, after ad- 
justing the radial-mode size by a pair of lenses (MA), to 
a first reflecting spatial light modulator (SLMl) which 
generates the desired 0AM qutrit state. The hologram 
kinoform displayed on the SLMl for each 0AM state to 
be generated, in the first-order diffraction, is shown in 
Fig.[l] After SLMl, a half wave-plate (HWP) and a quar- 
ter wave-plate (QWP) are used to write the polarization 
qubit in the photon. Hence we are able to generate any 
hybrid qusix that is a product of a qutrit and a qubit. 

The qusix-carrying photon is then sent to the detection 
stage. This stage is composed of a polarization analysis 
set (HWP, QWP and a PBS) and a second spatial Hght 
modulator (SLM2) for converting in diffraction the 0AM 
state to be detected back into a Gaussian mode. The 
photon is finally coupled to a single mode fiber, to filter 
only this Gaussian mode, connected to a single-photon 
counter module (SPCM). To eliminate the Gouy phase- 
shift effects between different 0AM eigenstates occur- 
ring in free propagation, an imaging system (not shown 
in the figure) has been included to image the screen of 
SLMl onto the SLM2. All waveplates and SLMs were 
computer-controlled so as to allow for a fully automatic 
generation and measurement procedure. With this setup, 
it is possible to perform a projective measurement upon 
every possible separable state of polarization and 0AM. 

As a first test, we verified the MUBs properties by 
generating each qusix \tpi) among the 18 states of the 
MUBs and then projecting it onto all the 18 states \ipj). 
Figure [2] shows the resulting measured probability dis- 
tribution Pij — I {ipj\ipi) p, compared to the theoretical 
one Pj'-. For a quantitative comparison, we used the sim- 

ilarity parameter S — ^ '' p ^^ 'p/ , which is a natural 
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Figure 2: Experimental analysis of hybrid qusix pho- 
tonic states: Probability distribution resulting from all 
18 X 18 projections of each state within the three MUBs over 
all the others, comparing theoretical and experimental val- 
ues. According to theoretical predictions, we expect that the 
18 X 18 matrix can be divided into nine 6x6 blocks A^ , 
where the two indices m,n £ {1,11,111} label generation and 
detection bases, respectively. Blocks that correspond to pro- 
jection of one basis over itself {m — n) should be diagonal, i.e., 
(yl™)ij = 5ij. Other blocks, whose values represent the over- 
lap between states belonging to two different bases, should be 
flat, i.e., (A™),, J = 1/6, for m ^ n. 



generalization of the fidelity used to compare two wave- 
functions, finding S = (99.19 ± 0.04)%. As a second 
check of the quality of our hybrid qusix states, we re- 
constructed the density matrix of all the 18 states by 
quantum state tomography. Since we lack a complete 
set of MUBs in dimension six, we performed measure- 
ments in all possible product states obtained combining 
the three MUBs of the polarization space tt and the four 
MUBs in the 0AM space O, for a total of 72 projections. 
In Table |T] the resulting experimental fidelities of the 18 
MUBs states are reported. The overall mean fidelity was 
F = (98.51 ± 0.04)%. Moreover, Figure [s] shows the re- 
constructed density matrices compared to the theoretical 
ones for three representative qusix states, one for each 
MUB considered here. 



PURE-OAM QUSIX ENCODING AND 
CHARACTERIZATION 

Our second experimental implementation of qusix pho- 
tonic states has been based on the 0AM space only. Al- 
though the hybrid approach may offer advantages for cer- 
tain specific tasks [41J, an encoding in 0AM is in princi- 
ple suitable of extension to arbitrary dimensionality and 



Basis 



State 



1^)1 + 1) 
\H) |0) 

1^)1-1) 

|V')|+1) 

\V) |0) 

|y)|-l) 



II 



III 



\A}\a,} 
\A)\a2) 
\A}\as) 
|D)|ai) 
|D)|q2) 

\L)\Pi) 

\L)\P2) 
\L)\P3) 

\R) Wi) 
\R) \P2) 

\R) m 



Average Fidelity 



Fidelity 



0.986 ± 
0.982 ± 
0.986 ± 
0.988 ± 
0.980 ± 
0.983 ± 



0.002 
0.002 
0.002 
0.002 
0.002 
0.002 



0.989 ± 
0.981 ± 
0.986 ± 
0.989 ± 
0.982 ± 
0.980 ± 
0.981 ± 
0.981 ± 
0.979 ± 
0.977 ± 
0.972 ± 
0.970 ± 



0.001 
0.002 
0.002 
0.001 
0.002 
0.002 
0.002 
0.002 
0.002 
0.002 
0.002 
0.002 



0.9851 ± 0.0004 



Table I: Experimental fidelities measured for all 18 qusix hy- 
brid states that characterize the three chosen MUBs. 



enables the generation of any kind of state, including the 
entangled ones which, for hybrid encoding, would need a 
more complex experimental setup. To define a 6D Hilbert 
space, we adopted the following 0AM eigenstates as log- 
ical basis: 



/ = {|-3),|-2),|-1),|1),|2),|3)}. 



(8) 



The three MUBs were still defined starting from the ten- 
sor products of a 2D and a 3D spaces, as given in Eq. 
(l5|. More details about the resulting states of the three 
bases /, //, /// are given in Methods. 

The experimental setup used for generating and testing 
the states of the MUBs is the same as in the hybrid qusix 
case (see Fig. ??), but with the polarization optics set so 
as to keep a fixed polarization everywhere. The kinoform 
generation was based on the method described in Sec. III. 
Figure |^a shows the intensity and phase profile of the 
18 0AM modes which form the three MUBs. In Figure 
|4]-b, the theoretical and experimental probability distri- 
butions for all combinations of state preparation and de- 
tection are reported. The similarity between the two dis- 
tributions is 5" = (99.06 ± 0.04)% while the mean fidelity 
over the 18 states is i^ = (98.78±0.08)%. Comparing this 
result with the hybrid case, in which only 0AM states in 
dimension 3 were generated, we find that the fidelity of 
the 0AM generation does not decrease rapidly with the 
dimensions. Hence, the holograhic method used in this 
work promises to be suitable for the high-fidelity genera- 
tion of 0AM photonic qudits with very large dimension 
d. 
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Figure 3: Quantum tomography of hybrid qusix photonic states: Density matrices associated to states of each of the 
three MUBs have been fully reconstructed by projections over all the 72 states obtained by direct product of the three MUBs of 
the 2D polarization space vr and the four ones of the 3D 0AM subspace C Here we show one state for each MUB. Experimental 
and theoretical matrices are reported for comparison. 
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Figure 4: Experimental analysis of pure OAM qusix: a) Graphical representation of all 18 states of the three selected 
MUBs, in the case of pure OAM 6D encoding. The precise definition of these states is given in Methods. For each state, both 
the intensity and phase patterns are shown, b) Theoretical and experimental probability distributions for an experiment in 
which all the 18 x 18 combinations of generated/detected states belonging to the three MUBs are tested. 



DISCUSSION 

In summary, we have reported the experimental imple- 
mentation of a non-extendable set of three MUBs for a 
photonic quantum system of dimension six by two differ- 



ent approaches. In the first, the qusix states have been 
implemented via a hybrid scheme based on polarization- 
orbital angular momentum encoding. All the 18 states 
belonging to the MUBs are in this case separable states of 
two different degrees of freedom. The demonstration of 



MUBs with high fidelity and unbiasedness has required 
the development of a new method for determining the 
kinoform to be visualized in the spatial light modula- 
tor. The second demonstrated approach was based on 
a quantum encoding in the photon 0AM space only, at 
a fixed polarization. The generation of a set of MUBs 
with high fidelity was again verified and this method is 
suitable for a convenient extension to higher dimension- 
ality. The techniques we have demonstrated here can find 
application in fundamental tests of quantum mechanics, 
quantitative wave-particle duality, quantum key distribu- 
tion, and tests of quantum complementarity and logical 
indeterminacy. 



tation: 



Methods 



In dimension d = 2, the eigenstates of the three Pauli 
operators provide a complete set of MUBs, which can be 
represented by the columns of the following three matri- 
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The 18 columns of these three matrices give the coeffi- 
cients of the logical basis superpositions defining the 18 
0AM states shown in Fig. pla). 
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In d = 3, there exist four MUBs. We represent them 
here as the columns of the following four matrices: 



This work was supported by project HYTEQ-FIRB, 
the Finanziamento Ateneo 2009 of Sapienza Universita 
di Roma, the European Research Council (ERG) Start- 
ing Grant 3D-QUEST, and by the Future and Emerg- 
ing Technologies (FET) program of the European Gom- 
mission, under FET-Open Grant No. 255914, PHOR- 
BITEGH. 



o.= 




0-2 = 



,©4 = 



73 



^/3 



1 1 


1 


1 UJ 


L^2 


1 ^2 


UJ 


1 1 


1 


aj2 UJ 


1 


^2 1 


UJ 



,(10) 



where uj = exp (i27r/3). 

In d = 6, we may construct three MUBs by a direct 
product of the tti, 7r2, tts bases and the corresponding 
first three bases Oi, O2, O^: 



I^TTl(E)Ol, II = TT2(E>02, /// = 713 (8)03. (11) 



These three 6D bases have the following matrix represen- 
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